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The transitions between the different vortex states of thin 
mesoscopic superconducting disks and rings are studied using 
the non-linear Ginzburg-Landau functional. They are saddle 
points of the free energy representing the energy barrier which 
has to be overcome for transition between the different vortex 
states. In small superconducting disks and rings the saddle 
point state between two giant vortex states, and in larger 
systems the saddle point state between a multivortex state 
and a giant vortex state and between two multivortex states is 
obtained. The shape and the height of the nucleation barrier 
is investigated for different disk and ring configurations. 

74.60.Ec, 74.20.De, 73.23.-b 



I. INTRODUCTION 

The study of superconducting samples with sizes com- 
parable to the penetration depth (A) and the coherence 
length (£) became possible due to recent progress in 
nanofabrication technologies. This evolution resulted in 
an increase of interest in the investigation of flux pene- 
tration and flux expulsion in such mesoscopic samples in 
order to explain the hysteresis behavior and the different 
phase transitions in thin superconducting samples [p]-^5| . 

It is well known that for type- II (ft = A/£ > l/v2) su- 
perconductors the triangular Abrikosov vortex lattice is 
energetically favored in the magnetic field range H c i < 
H < H c2 where n is the Ginzburg-Landau (GL) param- 
eter, and H c i and H C 2 are the first and second crit- 
ical fields of a Type-II superconductor. Since the ef- 
fective London penetration depth A = X 2 /d increases 
considerably in thin films one expects the appearance 
of the Abrikosov multivortex state even in thin Type-I 
(ft < l/v^) superconductors when the thickness d <C A. 
But, in small confined systems there is a competition 
between the boundary of the sample, which tries to im- 
pose the symmetry of the sample boundary on the vor- 
tex configuration, and this triangular Abrikosov state. 
As a consequence, the effective GL parameter k* = A/£ 
is no longer the only controlling parameter which deter- 
mines the shape of the vortex configuration in thin meso- 
copic superconducting samples ||. Previous theoretical 
and experimental studies of superconducting disks and 
rings [0-^6) found that, as function of the applied field, 
there are transitions between circular symmetric vortex 
states (called giant vortex states) with different vorticity 



L. Experimentally it was found that the magnetic field 
at which the transition L — > L + 1 occurs does not nec- 
essarily coincides with the magnetic field H tr where the 
vorticity of the ground state changes from L to L + 1, i.e. 
it is possible to drive the system in a metastable state. 
This is typical for first order phase transitions. For in- 
creasing applied field, the state with vorticity L remains 
stable up to the penetration field H p > H tr and transits 
then to the superconducting state with vorticity L + 1. 
For decreasing applied field, the state with vorticity L + l 
remains stable down to the expulsion field H e < H tr be- 
fore going to the state with vorticity L. This hysteresis 
effect is a consequence of the presence of an energy barrier 
between the states with vorticity L and L+l. The latter 
corresponds to different minima of the free energy in con- 
figurational space and the lowest barrier between those 
two minima is a saddle point. This barrier arises from the 
fact that the superconducting current around a vortex is 
in the opposite direction to the screening currents at the 
surface of the sample []l7|]. This Bean-Livingston model 
has been refined to different sample geometries [[l8] p3| . 
The time of flux penetration and expulsion is determined 
by the height of the energy barrier. 

The experimental consequences of the existence of 
these metastable states are: (i) hysteretic behavior [|J, 
(ii) paramagnetic Meissner effect p4|~p9||, (iii) fractional 
flux penetration 30 1 , and (iv) negative flux entrance Q , 
i.e. a decrease of the flux penetration through the super- 
conducting disk with increasing vorticity and increasing 
magnetic field. 

Schweigert and Peeters (3^] studied flux penetration 
and expulsion in thin superconducting disks and pre- 
sented an approach to find the saddle point states. They 
calculated the height of the free energy barriers which 
separate the stable states with different vorticity L. We 
will extend their approach and present a systematic study 
of flux penetration and expulsion in thin superconducting 
disks and disks with a hole in the center, i.e. mesoscopic 
ring structures. 

Bezryadin et al ]32j used the nonlinear GL equation to 
study the phase diagram of a thin-wire loop and a thin 
film with a circular hole in the limit k* 3> 1. They per- 
formed a stability analysis of the giant vortex state with 
vorticity L by allowing only the admixture of the L + l 
vortex state. A more rigorous stability analysis was per- 
formed by Horane et al |33j who studied the saddle points 
between two vortex states of a one dimensional wire of 
zero width. They allowed for more possible non-uniform 
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perturbations which may make the vortex state unstable. 
They found that the transition between two angular mo- 
mentum states occurs through a saddle point which has a 
zero in the order parameter at some point along the ring. 
Such a zero creates a phase slip center, allowing the phase 
winding required for the transition. Our systems have a 
non zero radial width and consequently such a scenario is 
not possible because the order parameter is not allowed 
to be zero along a radial line. In fact it was found in 
Ref. 1 31 1 that for a disk geometry, the saddle point for 
flux penetration corresponds to a state with suppressed 
superconductivity at the disk edge which acts as a nu- 
cleus for the following vortex creation. In the present 
paper we will find that for rings with a finite width this 
picture has to be modified because of the presence of two 
boundaries, i.e. two edges. 

Recently, Palacios [£6| calculated saddle point states 
and the energy barriers responsible for the metastabili- 
ties of superconducting mesoscopic disks using the low- 
est Landau level approximation. The central idea of 
his method was to find generic stationary solutions of 
the Ginzburg-Landau functional and to project the or- 
der parameter onto smaller subspaces spanned by a fi- 
nite number I of eigenfunctions, L2, where 
< Li < L2 < ... < Li. Palacios restricted himself 
to Z < 3 and therefore his approach is a special case of 
the one of Ref. |3l| where no such restriction on I was 
imposed and where also different radial states were in- 
cluded. Yampolskii and Peeters [Q investigated the in- 
fluence of the boundary condition (surface enhancement) 
on the superconducting states and the energy barriers 
between those vortex states. They also restricted their 
calculations to I < 3. 

Recently, Akkermans et al [^5| studied the behavior 
of metastable vortex states in infinite superconducting 
cylinders for k 3> 1, i.e. the London regime. They 
considered the situation where the vortices are symmet- 
rically distributed along a closed ring and they found 
structural phase transitions of vortex patterns between 
the metastable states. The key concept was the intro- 
duction of a special curve T, which embodies the main 
geometric features of a vortex configuration. This curve 
appears mathematically as a limit cycle of the system of 
currents generated by the vortex pattern and separates 
the paramagnetic and diamagnetic domains. 

The paper is organized as follows: In Sec. II we present 
the theoretical model and the calculation method to ob- 
tain the saddle points. In Sec. Ill we study thin super- 
conducting disks and extend and supplement our previ- 
ous results |3^| . We make a distinction between small 
and large superconducting disks. In small disks only the 
giant vortex state appears, while in larger disks multi- 
vortices can nucleate and transitions between different 
multivortices are possible |p5| , Pq| . In Sec. IV we consider 
superconducting rings, where we make a distinction be- 
tween small and large rings. Our results are summarized 
in Sec. V. 



II. THEORETICAL FORMALISM 

In the present paper we consider very thin supercon- 
ducting disks with radius R and thickness d, and super- 
conducting rings with inner radius Ri and outer radius 
R . These mesoscopic superconducting systems are im- 
mersed in an insulating medium in the presence of a per- 
pendicular uniform magnetic field Hq. To solve this prob- 
lem, we follow the numerical approach of Schweigert and 
Peeters For very thin disks and rings, i.e. Wd <C A 2 , 
with W = R the radius of the disk or W = R — Ri the 
width of the ring, the demagnetization effects can be ne- 
glected and the Ginzburg-Landau functional can be writ- 
ten as 



G — G r , 



(1) 



where G, G n are the free energies of the superconducting 
and the normal states, \& is the complex order parameter, 
a and (3 are the GL coefficients which depend on the 
sample temperature. L is the kinetic energy operator for 
Cooper-pairs of charge e* = 2e and mass m* — 2m, i.e. 



L = {-iti$ - e*^/c)/2m* 



(2) 



where A = 1? is the vector potential of the uni- 
form magnetic field Hq written in cylindrical coordinates 
p and <t>. 

By expanding the order parameter \P = Cifi in 
the orthonormal eigenfunctions of the kinetic energy op- 
erator L(fi — eiipi l^-Q, the difference between the super- 
conducting and the normal state Gibbs free energy can 
be written in terms of complex variables as 



F = G-G n 



(a + e,) QC* + ^A%C*C*C k C u (3) 



where the matrix elements AV, 



J dJ^tpiPjipkPi are cal- 
culated numerically. The boundary condition for these 
ifi , corresponding to zero current density in the insulator 
media, is 



A 



fi = 



(4) 



These eigenenergies q and the eigenfunctions ifi depend 
on the sample geometry. For thin axial symmetric sam- 
ples the eigenfunctions have the form ipj—t n! i)(p,(f) = 
exp(il(j))f n (p), where I is the angular momentum and the 
index n counts different states with the same /. Thus, 
the order parameter ^ can be written as 



l<Pn,l 



(5) 



We do not restrict ourselves to the lowest landau level ap- 
proximation (i.e. n = 1) and expand the order parameter 
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over all eigenfunctions with energy q < e* , where the cut- 
ting parameter is chosen such that increasing it does 
not influence the results. The typical number of complex 
components used are in the range N = 30 — 50. Thus 
the superconducting state is mapped into a 2D cluster of 
N particles with coordinates (xi,yi) <-> (Re(Ci), lm(d)), 
whose energy is determined by the Hamiltonian (^) . The 
energy landscape in this 2N+1 dimensional space is stud- 
ied where the local minima and the saddle points between 
them will be determined together with the corresponding 
vortex states. 

To find the superconducting states and the saddle 
point states we use the technique described in Ref. [ pl[ . 
A particular state is given by its set of coefficients {Ci}. 
We calculate the free energy in the vicinity of this point 
SG = G{C n ) - G(C) where {C n } is the set of coefficients 
of a state very close to the initial one. This free energy is 
expanded to second order in the deviations S = C n — C, 



SG = F m S* m 



where 



F m = (a + c-) C m + pAtfCjCtG , 
+ 2pA% n C k C? , 

D-mn = (3A k i CkC[ , 



(6) 



(7a) 
(7b) 
(7c) 



and I mn is the unit matrix. Using normal coordinates 
fim = XkQm we can rewrite the quadratic form as 
SG = 2 (jk%k + Vkxt) ■ To find the eigenvalues r/k and 
the eigenvectors "fk we solve numerically the following 
equation: 



B + Re(D) Im(L>) 
Im(L>) B - Re(D) 



Re(Q fc ) 
lm{Q k ) 



Vk 



Re(Q k ) 
lm(Q k ) 



(8) 



Starting from a randomly choosen initial set of coeffi- 
cients, we calculate a nearby minimum of the free energy 
by moving in the direction of the negative free energy 
gradient — 7fe. The set of coefficients of this minimum 
determines then the ground state or a metastable state. 
Starting from the initial set of coefficients we can also 
calculate a nearby saddle point state by moving to a min- 
imum of the free energy in all directions, except the one 
which has the lowest eigenvalue. In this direction we 
move to a local maximum. Repeating this procedure for 
many randomly choosen initial sets of coefficients {Ci} 
for fixed magnetic field, we find the different possible 
superconducting states and saddle point states. To cal- 
culate the magnetic field dependence we start from a su- 
perconducting state at a certain field and we change the 
applied field by small increments. By moving into the 
direction of the nearest minimum or saddle point, the 
corresponding state will be found for the new magnetic 
field, provided that the field step is small enough. 



III. SUPERCONDUCTING DISKS 

In the present section we discuss superconducting 
disks. Although the system is circular symmetric, in 
general we are not allowed to assume that ^(p,(f)) — 
F (p) e lL ^ because of the non-linear term in the GL func- 
tional. Nevertheless, in small disks the confinement ef- 
fects are dominant and this imposes a circular symme- 
try on the superconducting condensate, which means 
that only the 'giant' vortex state, i.e. a circular sym- 
metric vortex state, appears. For larger disks and not 
too large magnetic fields, the confinement effects are no 
longer dominant and multivortices can nucleate in a cer- 
tain magnetic field range. For this reason we make a 
distinction between small and large disks. 



A. Small disks: giant vortex state 

We consider superconducting disks with radius R = 
2.0£. First we investigate the influence of the number 
of terms in the expansion of Eq. ([5J) on the energy of 
the minima and the saddle points. For the approach 
that n = 1 (i.e. lowest Landau level) and if only one 
I is taken into account for each state, i.e. \& = Ciipi, 
we find three different states for three different values of 
i, I = 0, 1 and 2, in different magnetic field regions. In 
Fig. |l](a) the free energy F of these L-states, measured in 
the condensation energy F = a 2 itR 2 d/2(3, is shown by 
the dotted curves as a function of the applied magnetic 
field Ho. Next, we take into account two values of I 
and n — 1, i.e. \& = C^ipi^ + Ci 2 (pi 2 as was done in 
e.g. Ref. [g6|. With this approach we find L-states with 
L = 0, 1, 2 and because of the concomittant existence of 
two minima also saddle point states with I2) = (0, 1) 
and (I, 2) appear. These are the saddle point states for 
the transition between L-states with L = l\ and L = 
In Fig. 0(a) the L-states for this approach are given by 
the solid curves and the saddle point states by the dashed 
curves. The inset shows the transition between the L = I 
state and the L = 2 state in more detail. Notice that 
including one extra term in Eq. (^|) reduces appreciably 
the stability region of the different giant vortex states, 
i.e. its metastable region is strongly reduced. The L- 
states are only stable up to the point where its energy 
equals the saddle point states. 

Fig. |l|(b) shows the free energy as a function of the ap- 
plied magnetic field if we do not restrict ourselves to the 
lowest landau level and if we expand the order parame- 
ter over all eigenfunctions with energy < e* , where the 
cutting parameter e* is chosen such that increasing it 
does not influence the results. The solid and the dashed 
curves indicate respectively the stable L-states and the 
saddle point states. The transition between the L = I 
state and the L = 2 state is enlarged in the inset. Notice 
that the stability region of the L-states is further reduced. 
Allowing more basis functions in Eq. (||) does not have a 
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strong influence on the energy of the L-states, e.g. com- 
pare the dotted and the solid curves in Fig. 0(a), but 
it considerably decreases the energy of the saddle point 
between the L-states. In doing so, it reduces strongly 
the stability range of the metastable states, and conse- 
quently it reduces the size, i.e. the width in the magnetic 
field range, of the hysteresis effect p3|] . For example, 
we found (H tr /H c2 , H e /H c2l H p /H c2 ) a (1.0, 0.52, 1.25) , 
(1.0,0.715,1.245) and (1.0,0.73,1.24) for the L = <-> 
L = 1 transition when we include two, three and an 
arbitrary number of basis function in Eq. (|5j), respec- 
tively. Similarly, we found for the L = 1 <-> L = 2 tran- 
sition (H tr /H c2 ,H e /H c2 ,H p /H c2 ) a (1.715,1.52,1.81), 
(1.715,1.535,1.80), and (1.715,1.555,1.795) including 
two, three and an arbitrary number of basis functions, 
respectively. These results clearly show that one has to 
exert some caution to cutoff the expansion in Eq. (|^) 
when calculating the saddle point and thus the energy 
barriers. Notice that the expulsion field H e /H c2 is most 
strongly influenced by the number of terms in Eq. (||) . 

In Fig. H the transition barriers U, i.e. the energy dif- 
ference between the saddle point state and the nearby 
metastable states are plotted. We show the 'exact' nu- 
merical results (solid curves) and the results when in- 
cluding only two values of I with n — 1, i.e. \P = 
Ci x <fi x +Ci 2 (pi 2 (dashed curves). Notice that by approxi- 
mating the order parameter better, it substantially lowers 
the energy barriers, it increases the expulsion fields H e 
and lowers the penetration fields H p slightly. The en- 
ergy barrier is smaller for higher L — > L + 1 transitions 
which occur at larger magnetic fields. The inset shows 
the barrier between the L = 1 and the L = 2 state in 
more detail. 

The spatial distribution of the superconducting elec- 
tron density | V I'| 2 in the saddle point state corresponding 
to the transition from the L = 1 state to the L = 2 
states is depicted in Figs. ||(a-d) for the magnetic fields 
indicated by the open circles in the inset of Fig. |^, i.e. 
Hq/H c2 = 1.615, 1.665, 1.715 (i.e. the maximum of the 
barrier) and 1.765, respectively. In Ref. ]3l| ] similar re- 
sults were shown for the L = <-> L = 1 saddle point. 
High (low) Cooper-pair density is given by dark (light) 
regions. With increasing field, one vortex moves from 
the center to the outer region of the disk, and the state 
changes from L = 2 to L = 1. This is better illustrated 
by the contour plots of the phase of the order parameter 
which is shown in Figs. ^|(e-h) for the same configura- 
tions. Along a closed path, which lies near the edge of 
the superconductor, the phase difference Aip is always 
given by L times 2n, with L the vorticity or winding 
number. Light regions indicate phases ^ > and dark 
regions if < 2n. When encircling the superconductor 
near the boundary, we find that the phase difference Aip 
is equal to 2 x 2tt in Figs. |§(e-g) and Aip — 1 x 2ir in 
Fig. ||(h), which means vorticy L = 2 and 1, respectively. 
At the maximum of the barrier, i.e. when the energy of 
state L = 1 and L = 2 are identical, the saddle point 
transits from vorticity L = 2 to L = 1. At this point the 



Cooper-pair density is zero at the boundary of the disk 
which acts as a nucleation center for flux penetration and 
expulsion 

B. Large disks: multivortex states 

We consider now a larger superconducting disk with 
radius R = 4.0£ in which multivortex states can nucleate 
in certain magnetic field ranges [|6p5[| . Fig. [| shows the 
free energy as a function of the applied magnetic field 
Ho . The energy of the different L-states is given by solid 
curves when they are in the giant vortex state and by 
dashed curves when they are in the multivortex state and 
the saddle point states are given by the dotted curves. 
The open circles give the transition points between the 
multivortex state and the giant vortex state for fixed L. 
The inset shows the transition barrier U as a function of 
the applied field for the different L «-> L + 1 transitions. 
To distinguish qualitatively the giant vortex state from 
the multivortex state for fixed L we considered the value 
of the Cooper-pair density |\1/| 2 in the center of the disk. 
Fig. U shows |^| center which is zero for a giant vortex 
state and non zero in the multivortex states when there 
is no vortex in the center of the disk. For R = 4.0£ we 
find only multivortex states for L = 2, 3, 4 and 5 and the 
transition from the multivortex state to the giant vortex 
state occurs at H MG /H c2 = 0.52, 0.77, 0.875 and 0.935, 
respectively. Of course, for L = 1 there is no distinction 
between the giant and the multivortex state. For L = 5 
the spatial distribution of the superconducting electron 
density is given in the insets (a-d) of Fig. |^ at the mag- 
netic fields corresponding to the open circles in Fig. @, 
i.e. H /H c2 = 0.8, 0.85, 0.9 and 0.95, repectively. In the 
multivortex state the vortices move towards the center 
with increasing magnetic field and at the same time the 
vortices become wider, and therefore, the Cooper-pair 
density in the center decreases until the axial symmetry 
is recovered at the transition field Hmg = 0.935H c2 . 

Next, we will study the energy barriers U in more de- 
tail. For a superconducting disk with radius R = 4.0£ 
the energy barriers for the different L <-> L + 1 transi- 
tions are shown in the inset of Fig. |^. The height of the 
energy barrier for the L <-> L+l transition decreases with 
increasing L. The difference between penetration and ex- 
pulsion field decreases also with increasing L. Fig. |](a) 
shows the free energy of the L = 4 and the L = 5 
states in more detail (solid curves for the giant vortex 
state and dashed curves for the multivortex state) to- 
gether with the energy of the saddle point state between 
these states (dash-dotted curve). Fig. ^(b) gives the cor- 
responding energy barrier. The open circles correspond 
to the transition from multivortex to giant vortex state. 
The barrier height is clearly not influenced by the tran- 
sition from multivortex to giant vortex state, i.e. there 
are no jumps or discontinuities at the transition. The 
spatial distribution of the superconducting electron den- 
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sity l^j for this saddle point state is depicted in the 
insets of Figs, H(b) for the configurations indicated by 
the triangles, i.e. H /H c2 = 0.81, 0.885, 0.96 (the bar- 
rier maximum) and 1.035, respectively. Notice that also 
in the saddle point the transition between a multivortex 
state with L — 5 and a giant vortex state with L = 4 is 
clearly visible. 

Near the maximum of the barrier, the barrier height 
changes linearly with magnetic field. Therefore, we can 
approximate the energy barrier U near its maximum 



by 



U U max . ^ H Hmax 

Fq f h c2 



where the slope a is positive for H < H max and negative 
for H > Umax- In Fig. [?| the absolute value of the slope 
| a | is given as a function of L for H < H max by the 
closed circles and for H > H max by the open circles. 
The absolute value of the slope is different for the left 
and the right side of the maximum of the barrier. Notice 
that for L = and L = 1, |a| is larger for H < H max 
as compared to if > H max , while for L > 1 the reverse 
is true. For increasing L the slope decreases and the 
behavior could be fitted to 



(L)\ = T 



bL 



cL 



with a = 0.02586, b = -0.00300 and c = 1.40357 for H < 
H max , and a = 0.02322, b = -0.00217 and c = 1.26502 
for H > H max . These fitting curves are shown in Fig. |?j 
by the solid line for H < H max and by the dashed line 
for H > H max . In the inset of Fig. [?] the maximum 
of the barrier height U max is given by the symbols as a 
function of the vorticity L. The barrier height decreases 
for increasing vorticity and the behavior could be fitted 
to 



U„ 



<L) = 



bh 



1 



zsfh 



with a = 0.07229, b = -0.00791 and c = 0.48657, which 
is shown by the solid curve. 

For larger superconducting disks and higher values of 
L, different configurations of multivortices can occur with 
the same vorticity |2^,^7],|36| . Fig. || shows the free energy 
as a function of the applied field for the superconduct- 
ing states with vorticity L — 6 and L = 7. For both 
vorticities two configurations are possible; (i) L vortices 
on a ring and no vortex in the center (solid curve) and 
(ii) L — 1 vortices on a ring and 1 in the center (dashed 
curve). The vortex state is completely determined by 
the number of vortices in the center L center and the total 
number of vortices L. For this reason we characterized 
the states by the indices (L center ;L) in Fig. |8[ The in- 
sets (a-d) show the Cooper-pair density at H /H c2 = 0.6 
for the (0; 6) state, the (1; 6) state, the (0; 7) state, and 
the (1; 7) state, respectively. Notice further that for such 
large radius, there is no transition from a multivortex to a 



giant vortex state for these values of L. In Fig. ||the free 
energy of the saddle point state between the (1; 7) state 
and the (0; 6) state is given by the dash-dotted curve. 
This state describes the expulsion of one vortex when 
the system transits from the L = 7 to the L = 6 config- 
uration and is illustrated in Figs. ^ (a-c) where we show 
the spatial distribution of the superconducting electron 
density |\1/| 2 in the saddle point state at Hq/H C 2 = 0.5, 
0.6 and 0.7, respectively. To transit from L = 7 to L = 6, 
one vortex on the ring moves towards the outside of the 
disk and the vortex in the center takes the free place on 
the ring. High (low) Cooper-pair density is given by dark 
(light) regions. 



IV. SUPERCONDUCTING RINGS 

Now, we will consider superconducting disks with ra- 
dius R Q with a hole in the center with radius i?^. For the 
same reason as in Sec. Ill we make a distinction between 
small and large systems. 

A. Small rings: giant vortex state 

As an example, we consider a superconducting ring 
with radius R Q — 2.0£ and hole radius Ri = 1.0£. In 
Fig. |h] the free energy is shown as a function of the 
applied magnetic field for the different L-states (solid 
curves) together with the saddle point states (dashed 
curves). We find giant vortex states with L = 0, 1, 2, 3, 4. 
Comparing this result with the result for a disk with 
radius R — 2.0£, more L-states are possible and the 
superconducting/normal-transition moves to larger mag- 
netic fields H]. The inset shows the energy barrier U 
for the transitions between the different L-states as a 
function of the difference between the applied magnetic 
field Hq and the L — > L + 1 transition field Hl^l + i. For 
increasing L, the height of the energy barrier and the dif- 
ference between the penetration and the expulsion field 
decreases. The energy barrier near its maximum can be 
approximated by U/F a = U max /F + a(H - H max )/H c2 , 
and we determined the slope oil^l+i\ «o^i = —0.8 for 
H < H max and 0.9 for H > H max , a^ 2 = -0.6 for 
H < H max and 0.75 for H > H max , a 2 ^ 3 = -0.3 for 
H < H m ax and 0.42 for H > H max , and a 3 ^ 4 = -0.013 
for H < H m ax and 0.038 for H > H max . The slope de- 
creases again for increasing L and the absolute value of 
the slope for H < H ma x is smaller than for H > H max for 
every L, although the difference is relatively smaller than 
in the previous disk case where we found ao^i = —0.31 
for H < H max and 0.44 for H > H max , and a x ^ 2 = 
-0.06 for H < H max and 0.1 for H > H max . 

Next, we investigate the 2^3 saddle point. At 
Hq/H c2 — 2.01 (expulsion field) and 2.535 (penetra- 
tion field) the saddle point state equals the giant vor- 
tex states with L = 3 and L = 2, respectively. The 
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transition between these two giant vortex states is illus- 
trated in Figs. [ll](a-d) which show the spatial distribu- 
tion of the superconducting electron density |^| 2 corre- 
sponding with the open circles in the inset of Fig. |l0| at 
Hq/H C 2 — 2.1, 2.2, 2.315 (i.e. the barrier maximum) 
and 2.4, respectively. High (low) density is given by dark 
(light) regions. With increasing field one vortex moves 
from inside the ring, through the superconducting ma- 
terial, outside the ring. From Fig. |llj(c) one may infer 
that the Cooper-pair density is zero along a radial line 
and that the vortex is, in fact, a sort of line. That this 
is not the case can be seen from the left inset of Fig. |l2| 
which shows the Cooper-pair density |^| 2 along this ra- 
dial line for H /H C 2 = 2.315. The Cooper-pair density 
in the superconducting material is zero only at the cen- 
ter of the vortex which is situated at x m i n /^ « 1.5 and 
| '3/ 1 is very small otherwise, i.e. |*| 2 < 0.01. In Fig. |l| 
the position of the vortex, i.e. of x m i n , is shown as a 
function of the applied field. Over a narrow field region 
the vortex moves from the inner boundary towards the 
outer boundary. From H /H c2 = 2.01 to 2.25 the center 
of the vortex is still situated in the hole but the vortex 
already influences the superconducting state (see for ex- 
ample Figs. |ll](a,b)). From H /H c2 = 2.36 to 2.535 the 
center of the vortex lies outside the ring, but it has still an 
influence on the saddle point (see for example Fig. |ll](d)). 
In the region Hq /H C 2 = 2.25 — 2.36 the center of the vor- 
tex is situated inside the superconductor. This is also 
illustrated by the contourplot (right inset of Fig. |l2|) for 
the phase of the order parameter at Ho/H C 2 = 2.315, 
corresponding with the open circle in Fig. [l2| When en- 
circling the superconductor near the inner boundary of 
the ring, we find that the phase difference Aip is equal to 
2 x 2-7T which implies vorticity L — 2. When encircling the 
superconductor near the outer boundary, we find vortic- 
ity L — 3. If we choose a path around the vortex (located 
at Xmin), the phase changes with 27r and thus L = 1. At 
the transition field (Hq/H C 2 = 2.315) the center of the 
vortex of the saddle point is clearly not situated at the 
outer boundary as was the case for superconducting disks 
(see for example Figs. |](c,g), |6](b), pjh) and Ref. [pl|). 
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To illustrate this more clearly, Figs 
radial position of the vortex during 
tween the Meissner state and the L — 1 state, and be- 
tween the L = 1 state and the L — 2 state for a supercon- 
ducting ring with radius R Q = 2.0£ and for several values 
of the hole radius, i.e. R l /£_ = 0.0, 0.5, 1.0 and 1.5. The 
open circles indicate the ground state transition fields. 
Only for the case of the disk without a hole the cen- 
ter of the vortex at the saddle point occurs at the outer 
boundary of the disk for the magnetic field at which the 
ground state changes from L to L + 1. When the disk 
contains a hole in the center there are two boundaries 
and the center of the above vortex is now located be- 
tween those two boundaries. For a small hole with radius 
Ri = 0.5£ the position of the vortex can be approached by 
the arithmetic mean of the inner and the outer radius, i.e. 



x m in/£, ~ {Ro + Ri)/2, and for a larger hole with radius 
Ri = 1.5£ by the geometric mean y/R a Ri. The transition 
field increases and the magnetic field range, over which 
the transition occurs, decreases with increasing L. No- 
tice that the transition field for the L = 1 <-> 2 transition 
for Rj/£, = 0.5 is larger than the one for = 0.0 (see 
Fig. [13(b)), which agrees with Fig. 4 of Ref. pi]. 



B. Large rings: multivortex states 

First, we consider superconducting rings with radius 
R = 4.0£ and hole radius R l = 1.0£. In Fig. || the free 
energy is shown as a function of the applied magnetic 
field. The different L-states are given by solid curves 
for giant vortex states and dashed curves for multivortex 
states, while the saddle point states are given by the dot- 
ted curves. The open circles correspond to the transition 
between the multivortex state and the giant vortex state 
for fixed L. These transitions occur at Hmg/H C 2 = 0.93, 
1.035 and 1.14 for L = 4, 5 and 6, respectively. Notice 
that for such a small hole in the disk the maximum num- 
ber of L, i.e. L = 10, is the same as for the disk case 
without a hole (see Fig. ||). The spatial distribution of 
the superconducting electron density |\E'| 2 is depicted in 
the insets (a-c) of Fig. [lj for the multivortex state with 
L = 4 at H /H c2 = 0.8, L = 5 at H /H c2 = 0.9 and 
L = 6 at H /H c2 — 1.0, respectively. High (low) Cooper- 
pair density is given by dark (light) regions. Notice that 
there are always L — 1 vortices in the superconducting 
material and one vortex appears in the hole, i.e. in the 
center of the ring. 

The energy barriers for the transitions between the dif- 
ferent L-states are shown in Fig. [l^ as a function of the 
applied magnetic field. By comparing this with the en- 
ergy barriers for a disk with no hole, we see that the 
barrier heights and the transition fields are strongly dif- 
ferent (see the inset of Fig. ||). Therefore we show in 
the insets of Fig. [l^ the maximum height of the energy 
barrier U max and the L <-> L + 1 transition field H tr as 
a function of L for superconducting disks with no hole 
(squares) and with a hole of radius Ri = 1.0£ (circles), 
2.0£ (triangles) and 3.0£ (stars). In all cases the height 
of the energy barrier decreases and the transition fields 
increases with increasing L. By comparing the situation 
with no hole and with a small hole with Ri = 1.0£, we see 
that the barriers for L < 1 arc higher for the disk with 
a hole with Ri = 1.0£ than for Ri = 0.0£, while they are 
smaller when L > 1. Notice also that the value of the 
L — > L + 1 transition field is sensitive to the presence of 
the hole with radius Ri = 1.0£ for small L and insen- 
sitive for larger L. The reason is that for small L > 
such a central hole has always one vortex localised inside 
which favors certain vortex configurations above others, 
while for larger L in both cases only giant vortices appear 
with sizes larger than the hole size and the presence of 
the hole no longer matters. For larger holes the energy 
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barrier decreases more slowly, because the free energy of 
the different L states shows a more parabolical type of 
behavior as a function of the magnetic field. The transi- 
tion field has a much smaller dependence on the radius 
of the hole. Notice that the transition field for Ri = 3.0£ 
is linear as a function of L for small holes and L < 9. 
This is in good agreement with the results in the narrow 
ring limit, where the transition between states with dif- 
ferent vorticity L occurs when the enclosed flux <p equals 
(L + 1/2)0 O {§. 

Next, we investigate the saddle point states in these 
superconducting rings. We make a distinction between 
different kinds of saddle point states; i) between two gi- 
ant vortex states, ii) between a multivortex and a giant 
vortex state, iii) between two multivortex states with the 
same vorticity in the hole and different vorticity in the 
superconducting material, and iv) between two multivor- 
tex states with the same vorticity in the superconducting 
material but different vorticity in the hole. The first sad- 
dle point transition was already described for the case of 
small superconducting rings (see Figs. |ll| and |l2]) . Next, 
we study the saddle point state between a multivortex 
state with L = 5 and a giant vortex state with L = 4 for 
the previous considered ring with radius R — 4.0£ and 
hole radius Ri = 1.0£. Figs. [l6|(a-f) show the Cooper-pair 
density for these saddle point states at Hq/H C 2 = 0.83, 
0.88, 0.93, 0.965 (i.e. the barrier maximum), 1.03 and 
1.06, respectively. High (low) Cooper-pair density is 
given by dark (light) regions. For increasing field one 
vortex moves to the outer boundary, while the others 
move to the center of the ring where they create a giant 
vortex state. Remark that the giant vortex state is larger 
than the hole and therefore it is partially situated in the 
superconductor itself. 

To study saddle point transitions between different 
multivortex states we have to increase the radius of the 
ring to favour the multivortex states. Therefore, we 
consider a ring with radius R — 6.0£ and hole radius 
Ri — 2.0£. Fig. [17] shows the free energy of multivortex 
states with L = 8 and L = 9. In both cases 3 vortices are 
trapped in the hole. The lower insets show the spatial dis- 
tribution of the superconducting electron density |^| 2 at 
the transition field Hq/H C 2 = 0.695 for L = 8 and L = 9. 
It is clear that there are only 5 and 6 vortices in the su- 
perconducting material, respectively. The free energy of 
these multivortex states is shown by solid curves, while 
the saddle point energy between these states is given by 
the dashed curve. Notice further, that there is no transi- 
tion from the multivortex states to the giant vortex states 
with L = 8 and 9 as long as these states are stable. The 
spatial distribution of the superconducting electron den- 
sity I "J 1 2 for this saddle point state is depicted in the 
upper insets at the magnetic fields H /H C 2 = 0.63, 0.695 
(the barrier maximum) and 0.76, respectively. For in- 
creasing field one vortex moves from the superconducting 
material to the outer boundary and hence the vorticity 
changes from L — 9 to L = 8. Notice that the vorticity 
of the interior boundary of the ring does not change. 



The fourth type of saddle point state to discuss is the 
L — ► L+l transition between two multivortex states with 
the same vorticity in the superconducting material but 
with a different vorticity in the hole. For i? /£ = 4 and 
R /£ = 6 we did not find such transitions regardless of 
the hole radius. This means that at least for these radii 
there is no transition between such states which describes 
the motion of one vortex from the hole through the su- 
perconducting material towards the outer insulator. 

Finally, we investigated the influence of the hole radius 
on the barrier for a fixed outer ring radius. Fig. [TJ(a) 
shows the maximum barrier height, i.e. the barrier height 
at the thermodynamic equilibrium L — > L + 1 transition, 
as a function of the hole radius Ri for a ring with radius 
R = 4.0£ for the transition between the Meissner state 
and the L = 1 state (solid curve) and for the transition 
between the L = 1 and the L — 2 state (dashed curve). 
For increasing hole radius, the barrier height of the first 
transition rapidly increases in the range Ri — 0.1£ to 
Ri = 1.5£ and decreases slowly afterwards. For a super- 
conducting disk with radius R — 4.0£ with a hole in 
the center with radius Ri = 1.5£ the maximum barrier 
height for the — ► 1 transition is twice as large as for a 
superconducting disk without a hole. The barrier height 
of the second transition first decreases, then rapidly in- 
creases in the range Ri = 0.6£ to Ri — 2.5£ and then 
slowly decreases again. In this case the maximum bar- 
rier height for a superconducting disk with a hole with 
radius Ri — 2.5£ is three times as large as for a super- 
conducting disk without a hole. Hence, changing the 
hole radius strongl y i nfluences the maximum height of 
the barrier. In Fig. |l^(b) we plot the characteristic mag- 
netic fields of the barrier as a function of the hole radius, 
i.e. the transition magnetic field H tr , the expulsion mag- 
netic field H e and the penetration magnetic field H p , for 
the — > 1 transition by solid curves and for the 1 — ► 2 
transition by the dashed curve. For the — > 1 transition 
the characteristic magnetic fields decrease with increas- 
ing hole radius. For the 1 -> 2 transition the characteris- 
tic magnetic fields first increase to a maximum and then 
decrease. This behaviour was described and explained in 
our previous paper (see e.g. Fig. 17 of Ref. H]). Notice 
that the position of the minimum in U max coincides with 
the position of the maximum in H tr . 

V. CONCLUSIONS 

We studied the saddle points for transitions between 
different vortex states for thin superconducting disks and 
rings. A distinction was made between small systems 
where the confinement effects dominate and larger sys- 
tems where multivortices can nucleate for certain mag- 
netic fields. At the entrance of the vortex into the su- 
perconducting material the superconducting density be- 
comes zero at a certain point at the edge of the disk or 
ring. Such a zero in the order parameter acts as a phase 
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slip center which allows the vorticity to increase with one 
unit. For the case of the ring the vortex may enter (or 
exit) the superconducting material from the inner bound- 
ary or from the outer boundary of the ring. 

We studied the transition between two giant vortex 
states with different vorticity L. One vortex moves 
through the superconducting material to the center of 
the disk or to the hole. During the transition the po- 
sition of this vortex in the superconductor can be de- 
termined very precisely, because the Cooper-pair density 
is exactly zero in the center of this vortex. The transi- 
tion between a multivortex state and a giant vortex state 
with different vorticity L is also described. One vortex 
leaves (enters) the superconductor while the other vor- 
tices move towards (away from) the center of the disk. 
For large enough disk/ring radii, we calculated the tran- 
sition between two multivortex states. We found such 
transitions between two multivortex states with different 
vorticity L in the superconductor but with the same vor- 
ticity in the center/hole. One vortex enters/leaves the 
superconductor while the other vortices rearrange them- 
selves. Transitions between different multivortex states 
with the same vorticity in the superconducting material 
but different vorticity in the hole were not found for the 
considered ring configurations, which means that transi- 
tions between such states do not occur in these particular 
cases. 

The maximum height of the energy barrier always de- 
creases for increasing L. Near the maximum, the barrier 
height decreases linearly and its slope at the left side 
(H < H max ) of the maximum is not equal to the slope at 
the right side (H > H max ). The barrier shape and height 
strongly depend on the radius of the hole in the center 
of the disk. 
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FIG. 1. The energy of the minima in the free energy F 
and the energy of the saddle points as a function of the 
applied magnetic field Ho for a superconducting disk with 
radius R = 2.0£. (a) When only one term is included, i.e. 
(n, I) — (1, 1) (dotted curve), when two Z-values are included, 
i.e. (n, h) and (n, fa) (solid curves), with the correspond- 
ing energy of the saddle point (dashed curves); (b) The gi- 
ant vortex energy (solid curves) and the saddle point energy 
(dashed curves) when an arbitrary large number of terms are 
included. The free energy is scaled with the condensation 
energy F = a 2 nR 2 d/2p. 

FIG. 2. The transition barrier U for transitions between 
different L-states for a superconducting disk with radius 
R — 2. Of when taking into account only two values of I 
(dashed curves) and for the numerical 'exact' result (solid 
curves). The inset shows the second barrier in more detail. 



FIG. 4. The free energy as a function of the applied mag- 
netic field Ho- The different L-states are given by solid curves 
when in the giant vortex states and by dashed curves when 
in the multivortex states, while the saddle point states are 
given by the dotted curves. The open circles correspond to 
the transitions between the multivortex state and the giant 
vortex state for fixed L. The inset shows the transition barrier 
U as a function of the applied field for the different L <-» L + 1 
transitions. 

FIG. 5. The Cooper-pair density |^| 2 in the center of the 
disk with radius R = 4.0f for L — 2, 3, 4 and 5. The in- 
sets (a-d) show the spatial distribution of the superconducting 
electron density for L — 5 at the magnetic fields correspond- 
ing to the open circles; i.e. Ho/H C 2 = 0.8, 0.85, 0.9 and 0.95, 
repectively. The transition from multivortex state to giant 
vortex state occurs at the transition field Hmg- 

FIG. 6. (a) The free energy of the L = 4 and the L — 5 
states (solid curves for giant vortex states and dashed curves 
for multivortex states) and the saddle point states between 
these states (dash-dotted curve) for a superconducting disk 
with radius R — 4.0£; and (b) the energy barrier correspond- 
ing with this transition. The open circles correspond with the 
transition from multivortex to giant vortex state for fixed L. 
The insets show the spatial distribution of the superconduct- 
ing electron density |^| 2 for the saddle point states indicated 
by triangles, i.e. at the magnetic fields Ho/ Hc2 = 0.81, 0.885, 
0.96 (the barrier maximum) and 1.035. It is the transition 
between a multivortex state with L = 5 and a giant vortex 
state with L — 4. 

FIG. 7. The absolute value of the slope |a| of the energy 
barrier for a superconducting disk with R — 4. Of as a function 
of L for H < H max (closed circles) and for H > H max (open 
circles). The inset shows the maximum barrier height as a 
function of the vorticity L. The solid and dashed curves are 
the results of a fit. 



FIG. 3. The spatial distribution of the superconducting 
electron density |^[ 2 (a-d) and the phase of the order pa- 
rameter (e-h) in the saddle point state corresponding to the 
transition from the L = 1 state to the L = 2 states in a su- 
perconducting disk with radius R = 2. Of for the magnetic 
fields indicated by the open circles in the inset of Fig. ^| 
Hq/H C 2 = 1.615 (a,e), 1.665 (b,f), 1.715 (i.e. the maximum 
of the barrier) (c,g) and 1.765 (d,h). High Cooper-pair den- 
sity is given by dark regions, low Cooper-pair density by light 
regions. Phases ip > are given by light regions and <p < 2n 
by dark regions. 



FIG. 8. The free energy F as a function of the applied mag- 
netic field Ho of the (0; 6) and the (0; 7) state (solid curves), 
the (1; 6) and the (1; 7) state (dashed curves), and the saddle 
point state (dash-dotted curve) between the (1; 7) and (0; 6) 
state for a superconducting disk with radius R = 6. Of . The 
insets show the Cooper-pair density of the (0; 6) state (a), the 
(1;6) state (b), the (0; 7) state (c), and the (1; 7) state (d) at 
the transition field Ho/H C 2 = 0.6. 

FIG. 9. The spatial distribution of the superconducting 
electron density |<E'| 2 for the transition between the L = 6 
state and the L = 7 state for a superconducting disk with ra- 
dius 7?/f = 6.0 at the applied magnetic fields Ho/H C 2 = 0.5 
(a), 0.6 (b) and 0.7 (c). 
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FIG. 10. The free energy for a superconducting ring with 
radius R = 2.0£ and hole radius Ri = 1.0£ as a function of 
the applied magnetic field for the different giant vortex states 
(solid curves) and for the saddle point states (dashed curves) . 
The inset shows the energy barrier U for the transitions be- 
tween different L-states as a function of the difference between 
the applied magnetic field Ho and the L — > L + 1 transition 
field Hl^l+i- 

FIG. 11. The spatial distribution of the superconducting 
electron density |^| 2 of the transition between the giant vor- 
tex states with L = 2 and L — 3 for a superconducting ring 
with R = 2.0£ and Ri = 1.0£ at H /H c2 = 2.1 (a), 2.2 (b), 
2.315 (c) and 2.4 (d). High density is given by dark regions 
and low density by light regions. 

FIG. 12. The radial position of the vortex in the saddle 
point for the 2^3 transition through the superconductor 
with radius R = 2.0£ and Ri — 1.0£. The left inset shows the 
Cooper-pair density along de a;-direction at Ho/H C 2 = 2.315, 
and the right inset is a contour plot of the phase of the order 
parameter at Ho/H C 2 = 2.315. 

FIG. 13. The radial position of the vortex for (a) the «-» 1 
and (b) the 1 <-> 2 saddle point transition as a function of the 
applied magnetic field for a superconducting ring with radius 
Ro = 2.0£ and Ri = 0.0, 0.5, 1.0 and 1.5f. The open circles 
indicate the transition fields. 

FIG. 14. The free energy for a superconducting ring with 
Ro = 4.0£ and Ri = 1.0£ as a function of the applied mag- 
netic field for the different L-states (solid curves for giant 
vortex states and dashed curves for multivortex states), and 
the saddle point states (dotted curves) . The open circles cor- 
respond to the transition between the multivortex state and 
the giant vortex state for fixed L. The inset shows the spatial 
distribution of the superconducting electron density |^| 2 for 
the multivortex state with L = 4 at Ho/H C 2 = 0.8 (a), L — 5 
at Ho/H c2 = 0.9 (b) and L = 6 at H /H c2 = 1.0 (c). High 
Cooper-pair density is given by dark regions, low Cooper-pair 
density by light regions. 

FIG. 15. The energy barrier U for the transitions be- 
tween the different L-states in a superconducting ring with 
Ro = 4.0£ and Ri — 1.0£ as a function of the applied mag- 
netic field. The insets show the maximum height of the energy 
barrier Umax and the transition field H tr as a function of L 
for rings with R a = 4.0£ and Ri = 0.0, 1.0, 2.0, and 3.0£. 



FIG. 16. The Cooper-pair density for the saddle point state 
transition between a multivortex state with L = 5 and a giant 
vortex state with L = 4 at H /H c2 = 0.83 (a), 0.88 (b), 0.93 
(c), 0.965 (d), 1.03 (e) and 1.06 (f). High Cooper-pair density 
is given by dark regions, low Cooper-pair density by light 
regions. 

FIG. 17. The free energy of the multivortex states with 
L = 8 and L = 9 (solid curves) and the saddle point state 
(dashed curves) between these multivortex states for a super- 
conducting ring with R = 6.0£ and Ri = 2.0£ as a function 
of the applied magnetic field. The lower insets show the spa- 
tial distribution of the superconducting electron density \ty\ 2 
at the transition field Hq/H C 2 = 0.695 for L = 8 and L = 9. 
The upper insets show the spatial distribution of the super- 
conducting electron density \ty\ 2 for the saddle point states 
indicated by the open circles, i.e. at Ho/H C 2 = 0.63 (a), 
0.695 (b) and 0.76 (c). High Cooper-pair density is given by 
dark regions, low Cooper-pair density by light regions. 

FIG. 18. (a) The maximum barrier height as a function of 
the hole radius Ri for a ring with radius R = 4.0£ for the 
transition between the Meissner state and the L = 1 state 
(solid curve) and the transition between the L = 1 state and 
the L — 2 state (dashed curve); and (b) the transition mag- 
netic field H tr , the expulsion magnetic field H e and the pen- 
etration magnetic field H p as a function of the hole radius Ri 
for the transition between the Meissner state and the L = 1 
state (solid curve) and the transition between the L = 1 state 
and the L = 2 state (dashed curve) . 
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